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^ ; ABSTRACT 

^ : 

^ , We investigate X-ray isophote twists created by triaxiahty differences between the 

I luminous stellar distributions and the dark halos in elliptical galaxies. For a typi- 

^ I cally oblate luminous galaxy embedded in a more prolate halo formed by dissipation- 
less collapse, the triaxiahty difference of AT ~ 0.7 leads to typical isophote twists of 

^ ■ {Aipx) — 16° lb 19° at 3 stellar effective radii. In a model which includes baryonic 

. dissipation the effect is smaller, with AT ~ 0.3 and (A^x) — 5° it 8°. Thus, ac- 

^SJ ■ curate measurements of X-ray isophote twists may be able to set constraints on the 

! interactions between baryons and dissipationless dark matter during galaxy formation. 

I The 30° X-ray isophote twist in the E4 galaxy NGC 720 cannot be reproduced by our 

I model, suggesting an intrinsic misalignment between the halo and the stars rather than 

' a projection effect. 
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1. INTRODUCTION 



Determining the intrinsic shapes of the luminosity and mass distributions of elliptical galaxies 
is an important unresolved problem. Both observational constraints from individual systems and 
predictions from theory demonstrate that the stellar and dark matter components of ellipticals are 
well described as triaxial ellipsoids. Parametric (e.g., Ryden 1992; Lambas, Maddox, & Loveday 
1992) and non-parametric (Fasano & Vio 1991; Tremblay & Merritt 1995; Ryden 1996) inversions of 
the ellipticity distribution of the stellar components rule out the hypotheses that elliptical galaxies 
are all oblate or all prolate. Because of the information lost in projection, inversions that include 
triaxiahty cannot determine a unique shape distribution, but the ellipticals seem to fall into two 
classes: bright {Mb < —20), triaxial, boxy galaxies and faint, roughly oblate, disky galaxies (Trem- 
blay &; Merritt 1996). This division is consistent with other structural, dynamical, and evolutionary 
evidence for two classes of ellipticals (e.g., Fasano 1991; Fasano & Vio 1991; Busarello, Longo, & 
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Feoli 1992; Capaccioli, Caon, & D'Onofrio 1992; Nieto, Poulain, &; Davoust 1994; J0rgensen & 
Pranx 1994; Kormendy & Bender 1996). 

Binney (1978, 1985; see also Contopoulos 1956) first noted that a triaxial galaxy with its 
net angular momentum vector aligned with its short axis would show a "kinematic misalignment" 
between the projected axes of the rotation and the light, and surveys of elliptical galaxy kinematics 
have found that such kinematic misalignments are common (Davies &; Birkinshaw 1988; Franx, 
Illingworth, &; Heckman 1989; Jedrzejewski k, Schechter 1989). However, a kinematic misalignment 
may also be caused by an intrinsic misalignment between the rotational and the short axes because 
the angular momentum vector of an equilibrium triaxial system can lie anywhere in the plane of 
the long and short axes (Heiligman &: Schwarzschild 1979). Alternatively, the galaxy may not be 
in equilibrium (due to a recent accretion event or to tidal forces from a satellite galaxy), or it may 
be a misidentified SO with a bar (see Merritt 1997a for a review of examples). A statistical analysis 
of kinematic misalignments in bright ellipticals (Franx, Illingworth, &; de Zeeuw 1991) found limits 
on the mean stellar triaxiality and the mean intrinsic misalignment angle of (T*) < 0.7 and (V'int) < 
45°, where the upper bounds for both parameters occurred only for solutions with two widely- 
separated peaks. Although degeneracies prevented strong constraints on T,,,, there were indications 
for a bimodal population of galaxies, with a large fraction of nearly-oblate, short-axis rotators, and 
a small fraction of nearly-prolate, long-axis rotators. Tenjes et al. (1993) and Statler (1994a,b, 
Statler &; Fry 1994) have devised methods which invert the projected streamlines of stellar orbits 
to determine the intrinsic shape of a galaxy, and their initial results also indicate the existence of 
a bimodal population (Statler, Dutta, k Bak 1997; Bak & Statler 1997). 

Theoretical models of halo formation generally produce flattened, prolate-triaxial halos (e.g., 
Bardeen et al. 1986; Frenk et al. 1988; White & Ostriker 1990; Aguilar & Merritt 1990; Katz 1991; 
Cannizzo & Hollister 1992; Cole & Lacey 1996). For example, AT-body simulations of a small sample 
of halos by Dubinski (1991, 1992) and Dubinski & Carlberg (1991) found that at 25 kpc they are 
very flat (mean short axis ratio of (cd) — 0.42 ± 0.06) and nearly-prolate ((Id) ~ 0.8 ± 0.2). The 
final halo shapes and angular momenta were very sensitive to the initial cosmological tidal field, 
and the angular momentum vectors were well-aligned with the short axes, with an average intrinsic 
misalignment of ('0int) — 26° ± 22°. Warren et al. (1992) found similar results ((cd) — 0.62 it 0.13 
and {T^\) 0.68 it 0.24 at 40 kpc), and found that the more massive halos were slightly flatter 
and more prolate. Again, the angular momentum vector was most commonly aligned with the 
short axis. We expect dissipative processes to produce baryonic galaxies with shapes that are 
drastically different from those of their parent halos (e.g., Katz & Gunn 1991; Udry 1993; Navarro 
k White 1994), but the concomitant effects on the shapes of the halos are unclear. Preliminary 
simulations suggest that halos become rounder and more oblate through interactions with baryons 
(Evrard, Summers, k Davis 1994; Dubinski 1994), probably due to the destruction and damping 
of halo box orbits. Dubinski (1994) found that the flatness and prolateness are reduced relative 
to a dissipationless model (cd ~ 0.6 versus 0.4, and ~ 0.4-0.5 versus 0.8, at 20 kpc). Late, 
major mergers of disk galaxies in compact groups may have produced some fraction of the present 
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population of ellipticals (e.g., Heyl, Hernquist, & Spergel 1994; Weil k, Hernquist 1996; Barnes & 
Hernquist 1996). Although the simulations have not carefully explored the effects of the mergers 
on the shapes of the dark halos, Weil & Hernquist (1996) suggest that the shapes of remnant 
halos are rather round and oblate-triaxial, having no correlation with the shapes of the remnant 
stellar distributions. Current theoretical and observational results thus suggest that the modestly 
triaxial-oblate luminous parts of elliptical galaxies are embedded in halos which are more prolate, 
so that there should be a strong gradient with radius in the triaxiality of the mass distribution. 

Direct measurements of the shapes of halos are far more difficult than measurements of the 
shapes of the stellar distributions. Careful examination of the kinematics of polar rings (Arnaboldi 
et al. 1993), gas disks (Lees 1991; Bertola et al. 1991; Franx, van Gorkom, k, de Zeeuw 1994; Plana 
& Boulesteix 1996; Bureau & Freeman 1997), and planetary nebulae (Hui et al. 1995; Mathieu, 
Dejonghe, & Hui 1996) have been used to constrain the shapes of both stellar and dark matter 
potentials. Unfortunately, such tracers are rarely found at large enough galactocentric radii to be 
helpful for probing the dark matter potential, and it is likely that there are systematic correlations 
between the shapes of halos and the presence of rings and disks. Gravitational lenses can also 
be used to directly measure the shapes of mass distributions, and the data appear to require 
very flat mass distributions. However, the number of lenses is small, and the quantitative effects 
of external tidal perturbations are not yet understood (King &: Browne 1996; Kochanek 1996; 
Keeton, Kochanek, &; Seljak 1997; Witt &; Mao 1997). The most promising candidate for directly 
measuring the intrinsic shapes of nearby halos is high-resolution mapping of the X-ray emission 
from hot gas in the halo potential. For example, Buotc & Canizares (1994, 1996b, 1997; hereafter 
BC94, BC96b, BC97) examined the radial profiles, position angles, and axis ratios of the X-ray 
isophotes of the bright, isolated E4 galaxy NGC 720, and found that the shape and the position 
angle of the potential cannot be produced by the stars. Thus, by using only geometric evidence and 
the assumption that the gas is in quasi-hydrostatic equilibrium, they showed that the stars must 
be embedded in a more massive, flatter, dark matter halo. For the SO galaxy NGC 1332, Buote & 
Canizares (1996a) reached a qualitatively similar conclusion. 

X-ray observations can also be used to constrain the triaxialities of halos by searching for X-ray 
isophote twists. The projection of a triaxial light distribution whose triaxiality varies with radius 
produces isophotes whose axis orientations vary with radius (see Mihalas &: Binney 1981). Isophote 
twists in the optical surface brightness are common (e.g., Williams &: Schwarzschild 1979; Leach 
1981), and can be used to constrain the shape of an individual galaxy (Fasano 1995), although 
care must be taken to rule out intrinsic axis twist caused by other factors (Fasano & Bonoli 1989; 
Nieto et al. 1992). X-ray emission from a triaxial system can show an analogous isophote twist — 
as pointed out by Binney (1978) for the case of a galaxy cluster — allowing the triaxiality of an 
individual galaxy's potential to be probed. In the simplest case, if a galaxy's luminous and dark 
matter distributions are intrinsically aligned and have constant but different triaxialities, then the 
projected axes of its X-ray isophotes will twist from small radii, where the potential is dominated 
by the stellar core, to large radii, where it is dominated by the dark halo. Similarly, as discussed 
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by BC96b, the radial triaxiahty gradients produced in halos by dissipative processes can result in 
X-ray isophote twists. Here we employ simple models to study the behavior of X-ray twists, and to 
determine how well they constrain the triaxiahty of halos. In §2.1, we use an analytic approximation 
to find the amplitudes of twists expected for different assumptions about the shapes of galaxies 
and halos. In §2.2, we numerically calculate and project the X-ray emission for a small sample 
of models to verify the predictions of the analytic model, and to examine the radial behavior of 
the twist more closely. In §3, we attempt to model the large twist observed in NGC 720 (BC94, 
BC96b), and we present our conclusions in §4. 



2. CHARACTERISTICS OF X-RAY ISOPHOTE TWISTS 

2.1. Analytic Projection of Triaxial Ellipsoids 

An ellipsoidal luminosity density uirn?) with constant axis ratios c/a and h/a depends only on 
the ellipsoidal coordinate = X"^ /a^ + /b'^ + Z'^ /c^ (hereafter, we will set a = 1 for simplicity). 
When it is projected in the spherical-polar coordinate direction {0,(p), with the observer's z-axis 
along the line of sight, its surface brightness can be expressed as a one-dimensional integral of the 
radial density profile, 

2 

I{x,y) = —= / i'{u'^ + m1)du, where 
V / ^0 

f = sm^ 6'(cos^ (f) + — p-^) + = J + + ^2/ ), 

cos^9, . 9 , cos^o!), sin^^ „ ^ . ,, 1,1 

A= — 5— (sm^^H -r-!-) + —p^, B = cos 6' sm 20(1 - t^-)^, 

, „ ,sin^ (/) 2 ,x 1 
and C = {—^ + cos^ (fy-^. (1) 

For ^ = 0°, the line of sight is along the short axis; for 9 = 90°, = 0°, it is along the long axis; and 
for 9 = 90°, (j) = 90°, it is along the intermediate axis. The surface brightness / is also ellipsoidal, 
with a constant axis ratio q and a constant minor axis position angle from the projection of the 
short axis. If the galaxy rotates about its short axis, if) is equivalent to the kinematic misalignment 
angle. The axis ratio and the position angle are 



,2 A + C-^B- + {A-Cf 



tan2V; = ^; 7^ = rr^r.^.o n.:^2 . _,2^^ I .■,^2n - (2) 



B _ Tcos9sm24> 
A-C ~ r(cos2 9 sin^ - cos^ 0) + sin^ 6'' 

Equations 1 and 2 were derived by Binney (1985). The axis ratio g(6, c, 9, 0) is a function of the 
projection angles and both axis ratios of the triaxial ellipsoid, but the position angle ^(T, 9, (f)) 
is a function of only the triaxiahty parameter T= (1 — 6^)/(l — c^) and the projection angles 
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(Pranx, Illingworth, &; de Zeeuw 1991). If we now consider two triaxial distributions whose three- 
dimensional axes are intrinsically aligned, but which have different triaxialities (Ti < T2), they are 
misaligned in projection by the angle A'i/', given by 

. „^ , (r2-Ti)sin2(9cos6'sin2<3!) 
sm2Aw = , where 

D1 = cos^ Q sin^ + [Tj (cos^ Q sin^ ^ - cos^ (^) + sin^ Qf . (3) 

If the two distributions have the same triaxialities {T\ = T2), then the misalignment angle is zero 
(A-^ = 0); and if one of the distributions is oblate (Ti = 0), then the misalignment angle is equal 
to the other distribution's kinematic misalignment angle (A^ = V'2)- 

These expressions for the misalignment angle are applicable to the case of a luminous galaxy 

embedded in a dark halo, where the two distributions have differing triaxialities, labeled by T^, 
and Tfj. The projected mass will not be aligned with the projected light, an effect which can be 
important in gravitational lens systems (Keeton, Kochanek, & Seljak 1997). To demonstrate the 
amplitude of the misalignment for this simple model, Figure 1 shows the misalignment angle as a 
function of the viewing angles, Aijj{9,(f)), for T* = 0.25 and = 0.75. Large misalignment angles 
(A^ ^ 45°) occur when the distribution is viewed near the major axis (sin2 ^ 0.05 • AT/r*) at 
intermediate latitudes (T* ^ sin^ 9 ^T^). 

Figure 2 shows the mean misalignment angle, (Aijj), and its dispersion, (ta^, averaged over all 
viewing angles {6, for all possible triaxialities (T*, Tj). The dispersion aAtp is comparable to the 
mean, and (A^) increases smoothly with AT, so that the largest mean misalignments occur when 
the system consists of a very oblate spheroid and a very prolate spheroid. Figure 2 also illustrates 
the misalignment angles expected for a naive model of the shapes of galaxies and their halos. If we 
assume a luminous galaxy with rough shape limits taken from model la of Franx, Illingworth, &; de 
Zeeuw (1991), and embed it in a dissipationless dark halo whose shape is taken from the models 
of Dubinski & Carlberg (1991), we find large typical misalignments, (Aip) ~ 18° it 20°. If we use 
the model for a dissipational halo from Dubinski (1994), the typical misalignment angle drops to 
(Alp) ~ 7° lb 11°. In both cases, misalignments much larger and much smaller than the mean are 
common. Thus, for both galaxy formation models, the projected mass will for some viewing angles 
show a substantial position angle twist from small radii (where the luminous matter is dominant) to 
large radii (where the dark matter is dominant). If the X-ray isophotes show the same asymptotic 
axes as the projected mass, (Aip) will also describe the typical asymptotic X-ray misalignment 
angles. A full numerical model is required to confirm this assumption and to determine where in 
radius the isophote twist occurs. 
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0° 20° 40° 60° 80' 

Fig. 1. — The misalignment angle Ai/; of two triaxial bodies, for the spherical-polar viewing angles 4>, 9. 
The contours are spaced in 5° increments from 5° to 85°, and the solution for A-^ = 45° is highlighted by a 
dot-dash contour. 




Fig. 2. — Contours of the mean misalignment angle (A^') (left) and its dispersion ctai/i (right) between two 
distributions of triaxiality T, and Tj. As an aid to the eye, rough limits from the literature for stellar and 
halo triaxialities (T*, Tj) are shown for a model of galaxy formation which includes dissipation (bottom) and 
one which does not (top); the squares represent the median values, and the dotted curves show the "one-tr" 
limits. 
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2.2. Numerical Projection of the X-Ray Emission 

We model the galaxy as a triaxial "r/-model" distribution of stars (Dehnen 1993; Tremaine 
et al. 1994) embedded in a triaxial softened isothermal dark matter halo potential, with the axes 
aligned. The total density, p{r) = p*{rnl) + pd{'nT^), consists of two ellipsoids of the form 

An ideal, completely-ionized, single-temperature, non-rotating gas in hydrostatic equilibrium has 
a density Pg{r) = pgo exp[— /xmp$(r)/fcBT'g] and an X-ray emissivity 

i^x(r) = a(Ax, Tg)pl{r) = ^^e-^i^rn,Hr)/k^T, 

(Cavalicrc & Fusco-Femiano 1978), where /urrip is the mean mass per particle, the emissivity constant 
a (when convolved with the detector's spectral response) is weakly dependent on the gas temper- 
ature Tg, and the potential <I>(r) is evaluated numerically. The asymptotic profile is vy^{r) oc r~^^, 
where (3 = A'KGpQs\pmp/2>k-BTg. We integrate the X-ray emission along the line of sight to find the 
surface brightness, Ix{x,y) = J^^ux{X,Y, Z)dz, and we find the axis ratios q{a) and the minor 
axis position angles ^(a) of the X-ray isophotes by finding the maximum of Jx at each semi-major 
axis a. 

As shown by Buote &; Canizares (1994, 1996a), for a gas with any arbitrary temperature profile, 
the three-dimensional shape of its X-ray emissivity, z^x, is the same as that of the three-dimensional 
potential, in which it lies (i.e., Tx = T<j,). Although the potential is rounder than the mass, and 
its axis ratios have a different radial dependence, the projected mass and the projected potential 
are aligned because they must satisfy the two-dimensional Poisson equation. For an isothermal 
gas, we find that the X-ray emission i/x is also roughly aligned with the projected mass (~ 1.5° 
typical difference) , and thus confirm our assumption that equation 3 and Figure 2 give an accurate 
asymptotic representation of X-ray isophote twists. Furthermore, by treating the X-ray axis ratios 
bx and cx as those of a perfect triaxial ellipsoid of X-ray emissivity, we can use the analytic 
expressions (1-2) to find the position angle Vx(a) at any radius to an accuracy of ~ 1°, and the 
projected axis ratio qx{a) to an accuracy of ~ 0.01. 

Figure 3 illustrates in detail the X-ray twists of a realistic galaxy with both stars and a dark 
halo, for several viewing angles. The stellar component is specified by a Hernquist (1990) model 
(r/ = 2 in equation 4). The relative masses and scale radii of the halo and the stars are set so that 
the rotation curve along the long axis is fiat. The scale radius ratio is s*/sd ~ 1.8. The mass ratio 
of the dark matter and stellar components within the stellar effective radius (.Reff — l-8s*) is thus 
~ 1, and within 8i?cfr is ^ 6. The surface density profiles of galactic X-ray emission can typically be 
fit by King-type /3-models (see Sarazin &; Bahcall 1977) with /? = 0.4-0.6 (Forman, Jones, & Tucker 
1985; Fabbiano 1989), and we set /3 = 0.4. The axis ratios of the stellar component are chosen to be 
K = 0.96 and c* = 0.75, corresponding to = 0.18. The axis ratios of the halo in the dissipationless 
scenario are = 0.58 and Cd = 0.42, corresponding to = 0.81; and in the dissipational scenario 
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they are = 0.85 and Cd = 0.62, corresponding to = 0.45. The minor axis position angle of 

the X-ray isophotes, ipx, matches roughly the position angle of the light, ip^:, at small radii, and 
that of the halo, Vdj large radii, with a gradual transition at intermediate radii. The radial 
behavior (e.g., the location and sharpness) of this transition is insensitive to the gas temperature 
parameter /3, and is somewhat sensitive to the viewing angles, the stellar exponent rj, the scale 
radii, and the mass ratio. The half-way point of the twist generally occurs at semi-major axis a ~ 
0.1-0.4i?eff for the dissipationless halo model, and at ^ 0.4-0. 6i?eff for the dissipational halo model. 
As predicted by the analytic description, there is a large scatter in the twist amplitude, depending 
on the viewing angles. For the dissipationless halo, the median asymptotic twist is A^x ~ 15°, 
and our small sample shows twists as small as 3° and as large as 45°. For the dissipational halo, 
the median twist is ~ 6°, and twists from 1° to 20° are seen. Note that although the halo begins to 
dominate the total mass only outside ReS, the X-ray isophote twist can occur much closer to the 
core. This is because the triaxiality of the X-ray emissivity Tx is not a linear function of the axis 
ratios (bx, cx), and in this case where the stellar potential is nearly oblate, the small amount of 
additional flattening that is produced by the halo along the intermediate axis is sufficient to alter 
the triaxiality (and thus the isophote position angle) considerably. With our numerical model, we 
reappraise the limits from Figure 2 and find that at 3i?eff) the mean misalignment angle for the 
dissipationless model is (A^x) — 16° ± 19°, and for the dissipational model is ~ 5° it 8°. 

The axis ratios of the X-ray isophotes are less strongly affected by triaxiality differences and 
parameter adjustments. For the dissipationless halo, qx ~ 0.85-0.9 at i?efr for most viewing angles, 
and for the dissipational halo, qx ~ 0.9-0.95. Figure 4 shows X-ray isophote contours in the 
central regions of the galaxy for both halo models and some different viewing angles. Note that 
the larger twists occur for rounder isophotes, whose position angles are more difficult to determine. 
Although it is not apparent from the figure, the position angle of the X-ray isophotes asymptotically 
approaches that of the projected halo mass. In many cases we will not be able to resolve the X- 
ray twist, but we can still recognize its presence and measure its amplitude by observing the 
misalignment between the stellar and the outer X-ray isophotes. 

As a further complication, motivated by observations of the galaxy NGC 720 (BC97) , we add 
a contribution from discrete X-ray sources whose luminosity distribution is directly proportional 

to the stellar mass distribution , and whose total emission is 20% of the gaseous X-ray emission 
within 5i2cfr- This discrete component significantly modifies the triaxiality of the X-ray emission, 
and causes the X-ray twist to remain relatively flat to a much larger radius (see Figure 3). The 
character of the twist also becomes much more sensitive to the galaxy model parameters, including 
the gas temperature /?, once we include a discrete component. 
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FlG. 3. — The minor-axis position angle of the X-ray isophotes of a model galaxy (measured relative to the 
projected kinematic axis), as a function of the semi-major axis, for a gridded sample of viewing angles. The 
solid curves correspond to the dissipationless model, and the dashed curves to the dissipational model. Some 
special cases are also displayed for the dissipationless model (dotted curves). The center plot shows a case 
with a discrete source component added. The center left plot shows cases with different stellar component 
masses, M*. The top left plot shows cases with different values for the stellar exponent, rj. The top center 
plot shows the analytic approximation for the twist angle. The stellar and halo misalignment angles (V'* and 
^d) are marked (see top center plot for labeling). 
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Fig. 4. — X-ray isophotes of two model galaxies, for different viewing angles. The contours are logarithmi- 
cally spaced, with the central levels omitted for clarity. The major and minor axes of the projected mass are 
marked for the halo (dotted axes) and for the light (dot-dash axes). The kinematic axis \l/k is marked by a 
vertical solid line. The (twisting) X-ray major axis is shown by a solid line. The top two cases correspond 
to a model with a halo from the dissipationless scenario, while the bottom cases are from the dissipational 
scenario. The stellar effective radius is ~ 5.4 kpc. 
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3. A KNOWN EXAMPLE: NGC 720 

Finally, we consider the particular case of NGC 720, an isolated E4 galaxy with bright X-ray 
emission (Forman, Jones, &; Tucker 1985), and no evidence that environmental effects distort its 
X-ray isophote shapes significantly (BC96b). According to Binney, Davies, k, Illingworth (1990), 
the stellar misalignment angle is consistent with zero. Buotc &: Canizarcs (BC94, BC96b, BC97) 
have carefully studied and modeled NGC 720 based on data from ROSATs PSPC and HRI, and 
ASCA. They found that the X-ray isophotes arc significantly elongated (gx ~ 0.75 for a ~ 60"-110", 
where 1" ~ OAhgo kpc), with some indication that the ellipticity decreases at large radii (gx ~ 
0.8-0.9 for a ^ 120"). There is a significant X-ray isophote twist between the the inner and the 
outer isophotes, where the X-ray major axis position angle ^'x is consistent with being aligned 
with the optical angle at small radii (\l'x — 142° for a ^ 60"), and twists to a different angle 
at large radii (^x — 114° for a ^ 90"). The emission comes primarily from hot isothermal gas, 
with a lesser component (~ 20%) from X-ray binaries (BC97). BC94 demonstrated that the shape 
alone of the X-ray emission necessitates the presence of an extended dark matter halo, since the 
potential produced by the stellar distribution is too round to produce the high ellipticity observed 
for the X-ray emission. BC96b suggested triaxiality of the halo as a possible culprit for the X-ray 
isophote twist, but they were unable to reproduce the twist using a simple triaxial model in which 
the intrinsic axis ratios vary with radius according to a power law. 

We reinvestigate the hypothesis that triaxiality is responsible for the observed X-ray isophote 
twist by searching for a reasonable galaxy model whose stellar and halo components produce a 
radially- varying axis ratio gx(ffl) and minor axis position angle V'x(o) consistent with those re- 
ported for NGC 720. Taking advantage of the fact that we can use the analytic equations (1-2) 
to find good approximations for qx and ^x, we try to find reasonable values for the parameters 
^ = (r*, c*, Td, C(j, 0, ^) which will give us the observed data S = {ip*,q*, V'x, 9x)- We use a Bayesian 
statistical technique to find the most probable solution, and examine first the case with a dissi- 
pationless halo. For our priors, we use Gaussian distributions of the parameters (c*, cx, ?x), 
where (c*) = 0.75 ± 0.2, (cx) = 0.79 ± 0.08, (T*) = 0.2 ± 0.2, and (Tx) = 0.7 ± 0.3. For the data, 
we set = 0° ± 10°, g* = 0.60 ± 0.05, = 27° ± 8°, and gx = 0.75 ± 0.07, and we assume that 
the kinematic axis is intrinsically aligned with the short axis. 

Figure 5 illustrates the solutions, showing the total posterior probability of the parameters 
given the data, P{CM) = ^j^iP{i\S)/^jP{i\S), where P{i\6) = P{i)P{6\i). The maximum like- 
lihood occurs at {9 ~ 67°, (f) ~ 42°, c* ~ 0.51, cx ^ 0.71, ~ 0.17, Tx c± 0.96), corresponding 
to ■0* — 2.3°, V'x — 21.5°, ~ 0.62, gx — 0.81. This solution is a marginal fit to the data given 
the errors in the axis ratios and position angles of the X-ray isophotes. The solution for the stellar 
shape is reasonable — although the optical ellipticity of NGC 720 varies considerably with radius, 
there is no significant isophote twist (Jedrzejewski, Davies, & Illingworth 1987; Capaccioli, Piotto, 
Sz Rampazzo 1988; Peletier et al. 1990; Nieto et al. 1991; Sparks et al. 1991), indicating either that 
the stellar distribution is very oblate, or that it is triaxial and we are observing it nearly "edge-on" . 
The X-ray emissivity shape corresponds to a flat, prolate halo with Cd — 0.34, Td ~ 1.0, which is 
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consistent with the models of BC94, where Cd ~ 0.2-0.5, assuming an edge-on viewing angle and 
either an oblate or a prolate halo. If we use priors on the halo shape from the dissipational collapse 
scenario, it is significantly more difficult to fit the data. 




Fig. 5. — The projected likelihood distribution for the model parameters to fit the data for NGC 720. The 
prior distributions are plotted, as are the posterior distributions which included the priors, and the posterior 
distributions which did not include the priors (excepting the prior on 0). The squares mark the maximum 
likelihood solution. 



We next use the maximum likelihood model to see if it can reproduce the shapes and orienta- 
tions of the X-ray isophotes in detail. Since the stellar effective radius is R^s 52" (Burstein et al. 
1987), we take the scale radius to be s* = 28'.'6. After BC94, we set the gas temperature to give a 
radial profile exponent of /3 = 0.5. Since BC94 find that inside 450", the ratio of the masses of the 
halo and the stars ^ 4 at 90% confidence, with no upper bound, we use a model with a mass ratio 
of ~ 0.7 inside Re and ~ 7 inside 450" 9Res). Prom the fits of BC94, we take a typical halo 
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scale radius to be = 59"6. Figure 6 shows the model X-ray isophotes, and Figure 7 compares the 
radial variation of the position angle and axis ratio of its projected X-ray emission to the data|^. 
Although the twist and the ellipticity at large radii are marginally compatible with the data, the 
twist transition is too gradual and too close to the core. Despite much experimentation with the 
model parameters, we could not significantly improve our fit to the transition. By adding a discrete 
source component to the X-ray emission, we can push the transition radius further outward, but it 
becomes far more difficult to fit the amplitude of the twist. The effects are similar if we compute 
the twist with the moment method (BC94), or if we allow for substantial rotation of the X-ray 
emitting gas in the central regions. We thus qualitatively confirm the tentative result of BC96b 
that it is very difficult for simple models with a radial triaxiality gradient to simultaneously fit both 
the high ellipticity and the large, abrupt twist of the observed X-ray isophotes. 



100 




-100 

-100 -50 50 100 

X (arcKOc) 



Fig. 6. — X-ray isophotes of our model of NGC 720. (Compare BC94 Fig. 5.) The contours are spaced 
logarithmically by a factor of two, with the central levels omitted for clarity. The dotted contour shows the 
stellar surface density isophote at the stellar effective radius. The major and minor axes of the projected 
mass are marked for the halo (dotted axes) and for the light (dot-dash axes). The kinematic axis is 
marked by a vertical solid line. The (twisting) X-ray major axis is shown by a solid line. 



^Note that our models measure ^'x and gx along isophotes, while the data measure them using a moment method 
(compare cm and eiso from BC94); these do not differ significantly unless there are large radial gradients. 
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FlG. 7. — The major axis position angle (left) and the axis ratio (right) of the X-ray isophotes as a function 

of the semi-major axis. Our best-fit model is shown, where the parameters are calculated on isophotes (solid 
curves) and using a moment method (dotted curves) , and the data are shown with their 90% confidence limits 
(for clarity, data points at the same radius are offset slightly in radius). The PSPC data are indicated by 
squares, and the ffl?/data by crosses. The position angle (\1/*) and the axis ratio (g*) of the stellar isophotes 
are marked by dashed lines (taken from Jedrzejewski, Davies & lUingworth 1987; Capaccioli, Piotto, & 
Rampazzo 1988; Peletier et al. 1990; Sparks et al. 1991), as is the approximate position angle of the outer 
X-ray isophotes. 



4. CONCLUSIONS 

Current observational evidence indicates that the stellar parts of elliptical galaxies are modestly 
flattened and close to oblate, while current theories for the formation of galaxies suggest that the 
halos are more flattened and more prolate. The extent of this shape difference depends on the 
interactions between baryons and dissipationless dark matter during galaxy formation, but all 
current models imply that the mass distribution of an early type galaxy increases in flatness and 
prolateness with radius. There is already some evidence from the flatness of X-ray isophotes (Buote 
& Canizares 1994, 1996a) and from gravitational lens models (King & Browne 1996; Kochanek 1996) 
that the typical mass distribution may be flatter than the himinosity distribution. A difference in 
triaxiality between luminous galaxies and their dark halos is detectable through misalignments 
of gravitational lenses with their stellar distributions, and is a natural explanation for the second 
shear axis that is necessary to fit lensing models, although intrinsic misalignments and external tidal 
shear sources are also viable explanations (Keeton, Kochanek, & Scljak 1997). Such a triaxiality 
difference also produces X-ray isophote twists (Binney 1978; Buote & Canizares 1996b). 
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We examined the behavior of X-ray isophote twists using several simple models of a luminous 

galaxy of constant triaxiality embedded in an intrinsically-aligned dark halo of a different triaxiality. 
A simple analytic approximation gives accurate estimates of the asymptotic position angles of the X- 
ray isophotes, but predicting the detailed radial behavior of the twist requires numerical simulations. 
For a reasonable model of a galaxy and halo, we find that the X-ray isophote position angle makes 
a gradual transition from the center of the galaxy, where it is aligned with the position angle of the 
optical isophotes, to the periphery, where it is aligned with the position angle of the projected halo 
mass. The "half-way" point of the twist occurs well inside the stellar effective radius (~ 0.1-0.6 
i?efr)) and is not detectable given present X-ray resolution limits, but the misalignment of the X- 
ray isophotes at large radii with the stellar isophotes is more easily observable. By examining the 
amplitude of the misalignment for a large population of galaxies, we can in principle distinguish 
between alternative models for the halo shapes. The very prolate halos predicted by dissipationless 
collapse simulations should produce mean misalignments of (Aipx) — 16° ± 19° at ~ 3i?efr, while 
the more oblate halos predicted by simulations which include baryonic dissipation should produce 
smaller misalignments (~ 5° it 8°). In practice, measurements of these twists will be difficult because 
the twist angles are small, the X-ray isophotes are fairly round {{qx) ~ 0.9), and the number of 
isolated, nearby X-ray galaxies is limited. Our results are not very sensitive to the parameters 
of the model, but the presence of a discrete source component can modify the detailed properties 
of the twist and should be included in any statistical test. Note that while we scaled our models 
to match the halo shapes predicted by simulations, any mass distribution with a strong radial 
triaxiality gradient will produce an X-ray isophote twist when observed from the proper angles. 
For example, the disruption of box orbits by a central stellar cusp can cause the central parts of 
an elliptical galaxy to become more oblate than the outer parts due to the longer time scales for 
scattering and phase mixing of the outer orbits (e.g. Merritt 1997b), leaving a twisting signature 
in the X-ray isophote position angles. 

We attempted to model the observed X-ray twist of the bright, isolated E4 galaxy NGC 720 

(BC94, BC96b), whose X-ray isophotes arc strongly flattened, and show a large, abrupt position 
angle twist of 30° at ~ 1-2 RcS- While a model halo from the dissipationless scenario fits 
the data better than does a halo from the dissipational scenario, it is in both scenarios difficult 
to reproduce both the twist behavior and the isophote flatness using reasonable parameters for 
the galaxy and halo. The simplest explanation for NGC 720 is that the halo and the stars are 
intrinsically misaligned. Such an intrinsic misalignment may arise naturally for a halo which forms 
by dissipationless collapse, or it may be caused by late-history major mergers. Whether caused by 
triaxiality or by intrinsic misalignment, the shapes and orientations of stellar and X-ray isophotes 
are important fossil clues to the formation history of galaxies. 



We thank David Buote and Christine Jones for their helpful comments. CSK is supported by 
NSF grant AST-9401722 and NASA ATP grant NAG5-4062. 
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